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Abstract 

Using an eikonal structure for the scattering amplitude, Block and Kaidalov[^] have derived factor- 
ization theorems for nucleon-nucleon, 7p and 77 scattering at high energies, using only some very general 
assumptions. We present here an analysis giving experimental confirmation for factorization of cross sec- 
tions, nuclear slope parameters B and p-values (ratio of real to imaginary portion of forward scattering 
amplitudes), showing that: 

• the three factorization theorems[Q| hold, 

• the additive quark model holds to ~ 1%, 

• and vector dominance holds to better than ~ 4%. 
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1 Introduction 



Assuming factorizable eikonals in impact parameter space b for nucleon-nucleon, 7p and 77 scattering pro- 
cesses whose opacities are equal, Block and KaidalovQ have proved three factorization theorems: 

where the cr's are the total cross sections for nucleon-nucleon, 7p and 77 scattering, 

Bnn{s) _ S-yp(s) 



2. 



3. 



where the B's are the nuclear slope parameters for elastic scattering. 



where the p's are the ratio of the real to imaginary portions of the forward scattering amplitudes, 

with each factorization theorem having its own proportionality constant. These theorems are exact, for all 
s (where y/s is the c.m.s. energy), and survive exponentiation of the eikonal||l|. 

Physically, the assumption of equal opacities, where the opacity is defined as the value of the eikonal at 
& = 0, is the same as demanding that the ratios of elastic to total cross sections are equal, i.e., 
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as the energy goes to infinity 

Factorization theorem ^ involving ratios of cross sections, is perhaps the best known. Factorization 
theorems ^ and |^ are less known, but turn out to be of primary importance. The purpose of this note is to 
present strong experimental evidence for all three factorization theorems, as well as evidence for the additive 
quark model and vector dominance. 



2 Eikonal Model 

In an eikonal model[^, a (complex) eikonal x(^5 •s) is defined such that a{b,s), the (complex) scattering 
amplitude in impact parameter space b, is given by 



a{b, s) = ^ (1 - e'^'""^) = ^ (1 - e-'^i(''^^)+*XH(^-)) . 
Using the optical theorem, the total cross section (Ttot(s) is given by 



atot(s) =2 / [l-e-Xi(''^^)cos(xH(6,s)) 



d^b, 



the elastic scattering cross section aci{s) is given by 



Colastic('S) 



and the inelastic cross section, di 



inelastic 



(s), is given by 



f inelastic (s) = O'tot(s) — Crelastic(s) — 



1 _ e^^^i^''^'') 



The ratio of the real to the imaginary part of the forward nuclear scattering amplitude, p, is given by 

Rc(i(ri - e-Xi )d2 61 

P(^) = -A— 4 

Im |i(/(l - e-'^i(''^")+*XR('''^)) d26| 

and the nuclear slope parameter B is given by 
/ b^a{b,s)d^b 



B = 



2 / a(&, s)d^b' 



(2) 
(3) 
(4) 
(5) 

(6) 
(7) 
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2.1 Even Eikonal 



A description of the forward proton-proton and proton-antiproton scattering amplitudes is required which is 
analytic, unitary, satisfies crossing symmetry and the Froissart bound. A convenient parameterization!^, ^ 
consistent with the above constraints and with the high-energy data can be constructed in a model where the 
asymptotic nucleon becomes a black disk as a reflection of particle (jet) production. The increase of the total 
cross section is the shadow of jet-production which is parameterized in parton language. The picture does 
not reproduce the lower energy data which is simply parameterized using Regge phenomenology. The even 
QCD-inspired eikonal Xcvcn for nucleon-nucleon scattering|^, ^ is given by the sum of three contributions, 
gluon-gluon, quark-gluon and quark-quark, which are individually factorizable into a product of a cross 
section a{s) times an impact parameter space distribution function W{b;fi), i.e.,: 

- Xgg(s,6)+Xqg(s,5)-fXqq(s,6) 

= « [f^gg{s)W{b ; ^gg) -f aqg{s)W{b ; VMqqMgg) + Crqqis)W{b ; /iqq)] , (8) 

where the impact parameter space distribution function is the convolution of a pair of dipole form factors: 

w{b■,^^)^^{^ibfK,{^ib). (9) 

It is normalized so that J W{b; ^)(Pb — 1. Hence, the cr's in eq. (||) have the dimensions of a cross section. 
The factor i is inserted in eq. (|^) since the high energy eikonal is largely imaginary (the p value for nucleon- 
nucleon scattering is rather small). 

The opacity of the eikonal, its value at 6 = 0, is given by 

O"" = [c^gg('S)/^gg + 0-qg(s)Aiqg + 0-qq(s)/^qq] , (10) 

a simple sum of the products of the appropriate cross sections a with the /^^'s, a result which we will utilize 
later. 

2.2 Odd Eikonal 

The odd eikonal, Xodd{b,s) = i(7oddW{b; fio^ii), accounts for the difference between pp and pp, and must 
vanish at high energies. A Regge behaved analytic odd eikonal can be parametrized as (see Eq. (5.5b) of 
Ref. 1) 

Xoddib, S) = - Codd^gg'^e"/^ W{b; flodd) (11) 

where /iodd is determined by experiment and the normalization constant Codd is to be fitted. 

2.3 Total Eikonal 

The data for both pp and pp are fitted using the total eikonal 

Xpp ~ XovGii i Xodd- (12) 



3 A Global Fit of Accelerator and Cosmic Ray Data 

Using an eikonal analysis in impact parameter space. Block et al. |3[ ^ have constructed a QCD-inspired 
parameterization of the forward proton-proton and proton-antiproton scattering amplitudes which fits all 
accelerator data[^ for <Ttot, nuclear slope parameter B and p, the ratio of the real-to-imaginary part of the 
forward scattering amplitude for both pp and pp collisions, using a procedure and the eikonal of eq. (|T^); 
see Fig.0 and Fig. I which are taken from ref. 1^ — in addition, the high energy cosmic ray cross sections 
of Fly's Eye and AGASA |^ experiments are also simultaneously fit|^. Because the parameterization is 
both unitary and analytic, its high energy predictions are effectively model-independent, if you require that 
the proton is asymptotically a black disk. A major difference between simultaneously fitting the cosmic ray 
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Figure 1: The fitted o-pp and crpp, in mb vs. y/s, in GcV, for the QCD-inspired fit of total cross section, B and p for both pp 
and pp. The accelerator data (squares are pp and circles are pp ) and the cosmic ray points (diamond, Fly's Eye and triangles, 
AGASA) have been fitted simultaneously. The cosmic ray data that are shown have been converted from cJJ^l^j^ to o-pp. 

and accelerator data and earlier results in which only accelerator data were used, is a lowering (by about a 
factor of w 2) of the error of the predictions for the high energy cross sections. In particular, the error in 
CTpp at y/s — 30 TeV is reduced to « 1.5 %, because of significant reductions in the errors estimated for the 
fit parameters (for a more complete explanation, see ref. 1^). 

The plot of cTpp vs. ^/s, including the cosmic ray data, is shown in Fig.[^, which was taken from ref. j^]. 
The overall agreement between the accelerator and the cosmic ray pp cross sections with the QCD-inspired 
fit, as shown in Fig.|^, is striking. 

In brief, the eikonal description provides an excellent description of the experimental data at high energy 
for both pp and pp scattering at high energies. 

4 Factorization 

We emphasize that the QCD-inspired parameterization of the pp and pp data ||^, ^, ^ allows us to calculate 
accurately the even eikonal of eq. (H) needed for: 

• the total cross section CTnn (from eq. (^)) used in the factorization theorem |l|, 

• the nuclear slope parameter B^n (from eq. (0)) used in the factorization theorem ^, 

• and the p- value pnn (from eq. (|^)) used in the factorization theorem H, 
since we must compare nn to 7p and 77 reactions. 

4.1 Theorems 

As shown in ref. the eikonals for 7p and 77 scattering that satisfy eq. (||)are given by 



and 



Kagg{s)W{b ; \j -^gg) 4- KCTqg(s)VF(& ; W-V^qqAigg + Kcrqq{s)W{b ; \ -Hqq) 



KVgg(s)VF(& ; i/Xgg) + KVqg(s)VK(6 ; ^V^^qq^gg + KVqq(s)W^(6 ; ^/^qq) 



(13) 



(14) 
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Figure 2: The fitted values for the nuclear slope parameters Bpp and Bpp, in (GeV/c) ^ vs. y/s, in GeV, for the QCD-inspired 
fit are shown in (a). In (b), the fitted values for ppp and ppp are shown. 



where we obtain x'*^ from x''™" from multiplying each a in ^^^^^ by k and each fi by y , and, in turn, 

we next obtain x^'^ from from multiplying each a in x''^ by k and each ^ by -y/^- The k in eq. ( p^ ) 

and eq. (^) is an energy-independent proportionality constant. We emphasize that the same k must be 

used in the gluon sector as in the quark sector for the ratio of (^7) to be process- independent (see eq. (j^). 

The functional forms of the impact parameter distributions are assumed to be the same for 773, 77 and nn 
reactions. It is clear from using eq. (^, eq. (p^), eq. ( p^ ) and then comparing to the opacity of eq. (pO|), that 
the three opacities are all the same, i.e., 



Hence, from ref. Q, we have the three factorization theorems 



,(s) _ a^p(s) 



(16) 



(17) 



Pnn(s) ^ Pfpjs) 

valid for s. It is easily inferred from ref. |^ that 



(18) 



T^^(s) = K^'had'^7p(s) = ('«^had)^^n„(s) (19) 
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B-l-yis) 

Puis) 



= P7p(s) 



= Pnn(s), 



(20) 
(21) 



where P^^^ is the probabihty that a photon transforms into a hadron, assumed to be independent of energy 
and K is a proportionahty constant, also independent of energy. The value of k, of course, is model-dependent. 
For the case of the additive quark model, = §■ 

We emphasize the importance of the result of eq. ( pl| ) that the p's are all equal, independent of the 
assumed value of k, i.e., the equality does not depend on the assumed model. 

4.2 Experimental Verification of Factorization using Compton Scattering 

The solid curve in Fig. ^ is pnn, plotted as a function of the c.m.s. energy ^s. According to eq. (pT]), this 




-0.2 :/ 



100 
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Figure 3: The solid curve is Pnn, the predicted ratio of the real to imaginary part of the forward scattering amplitude for the 
'elastic' reactions , 7 + p — > y + p scattering amplitude, where V is p, u> or ij> (using the factorization theorem of eq. (|2l[)). The 
dotted curve is ratio of the real to imaginary part of the forward scattering amplitude for Compton scattering , 7 + p ^ 7 + p, 
found from dispersion relations It has been slightly displaced from the solid curve for clarity in viewing. 

should be the same as pjp. No experimental data for the 'elastic scattering' reactions 7 + p V + p, where 
V is the vector meson p, uj or (/), are available for direct comparison. However, Damashek and Oilman |^ 
have calculated the p value for Compton scattering 7 + p ^ 7 + P using dispersion relations, i.e., the true 
elastic scattering reaction for photon-proton scattering. The dispersion relation calculation gives p^p if we 
assume that it is the same as that for the 'elastic scattering' reactions j + p ^ V + p. In this picture we 
expect that B-yp = Bp = B^^ — B^. We then compare the dispersion relation calculation, the dotted line in 
Fig. H, with our prediction for p^p from eq. ( pi] ) (pnn, the sohd Hne taken from ref. Q). The agreement is 
so close that the two curves had to be moved apart so that they may be viewed more clearly. It is clearly 
of importance to extend the energy region of the dispersion calculation. However, over the limited energy 
range available from the dispersion calculation, the prediction from eq. (|2^) of equal p- values is well verified 
experimentally. 



4.3 Quark Counting 

The additive quark model tells us from quark counting that k in eq. (^0|) is given by k = |. We can 
experimentally determine k by invoking from eq. (|2^) the relation B^p = Ki?nn {Bnn is computed using the 
parameters from ref. Q), and fitting n. In our picture, the 'elastic scattering' reactions 7 -I- p — > V + p, 
where V is the vector meson p, lu or 0, require that Bp = B^^ = B^{— S-yp). To determine the value of n 
in the relation B^p = kB^^, a fit was made to the available B^p data. In Fig. ^ we plot kB^u vs. the 
c.m.s. energy ^/s, using the best-fit value of k = 0.661, against the experimental values of B-yp. The fit gave 
K = 0.661 ± 0.008, with a total — 16.4 for 10 degrees of freedom. Inspection of Fig. ^ shows that the 
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experimental point of Bp at y/s = 5.2 GeV — which contributes 6.44 to the — clearly cannot lie on any 
smooth curve and thus can safely be ignored. Neglecting the contribution of this point gives a x^/d.f.==0.999, 
a very satisfactory result. We emphasize that the experimental 7p data thus 

• require k = 0.661 ± 0.008, a w 1% measurement in excellent agreement with the value of 2/3 that is 
obtained from the additive quark model. 

• clearly verify the nuclear slope factorization theorem of eq. ( pO[ ) over the available energy range spanned 
by the data. 




Figure 4: A fit of experimental data for the micJear slopes B, from the 'elastic scattering' reactions 7 + p — > V + p, where 
V is p, u> or (p, to the relation B-yp = KBnn, of eq. (|20|), where k = 0.661 it 0.008 

For additional evidence involving the equality of the nuclear slopes Bp, B^ and i?^ from differential elastic 
scatttering data see Figures (13,14,15) of ref. [||. 

4.4 Vector Dominance, using 7p Cross Sections 

Using K = I and eq. (p^, we write 

^7p(^) = ^^had'^n„(s), (22) 

where i^^ad probability that a photon will interact as a hadron. We will use the value P^^^ = 1/240. 

This value is approximately 4% greater than that derived from vector dominance, 1/249. Using (see Table 
XXXV, pag. 393 of ref. |o)) /2/47r = 2.2, fl/An = 23.6 and /|/47r = 18.4, we find Ey(47ra//2 ) = 1/249, 
where V = p,u}, 4>. The value we use of 1/240 is found by normalizing the total 77? cross section to the low 
energy data and is illustrated in Fig. |5[ where we plot the total cross section for 7 + p ^ hadrons from 
eq. (|2^) as a function of the c.m.s. energy -/s. The values for dnn have been deduced from the results of ref. 
Q , using the even eikonal from eq. (|8|) . The fit is exceptionally good, reproducing the rising cross section for 
7P, using the parameters fixed by nucleon-nucleon scattering. The fact that we use the value 1/240 rather 
than 1/249 (4% greater than the vector meson prediction) reflects the fact that P^^^^, the total probability 
that the photon is a hadron, should have a small contribution from the continuum, as well as from the vector 
mesons p, and uj. Thus, within the uncertainties of our calculation, the experimental data in the jp sector 

• are compatible with vector meson dominance. 

• agree with cross section factorization theorem of cq. (|l9|). 
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Figure 5: Th p tn tal cross section for 7p scattering. The solid curve is the predicted total cross section from the factorization 
relation of eq. ( |l9[ ) , cr-^p = | P]^^^^ ""nn , where P^^^ 



1/240. 



4.5 Experimental Verification of Factorization using 77 Scattering 



Using quark counting and the factorization theorem of eq. (|l^), we now write a.yy 

-Phad = 1/240. In Fig. I we plot our factorization prediction for a-yy{s) as a function of the c.m.s. energy 



l-Phad) (^nn ^hcrC 



factorization 

• Lep2-L3 

■ Lep2-0PAL 
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Figure 6: The predicted total cross section for 77 scattering from the factorization theorem of eq. (p^), cr-, 
1/240. The data sources are indicated in the legend. 



where P^^^ 



\ 3 had / 



■y/s and compare it to various sets of experimental data. It is clear that factorization, as expressed in 
eq. (p^, selects the preliminary L3 data (solid circles) rather than the preliminary OPAL results (solid 
squares) [|ll|. The Monte Carlo-averaged final results of L3, given by the open circles, agrees, within errors, 
with the revised OPAL data, with both new sets having a normalization of about 15-20 % higher than the 
factorization prediction given by the solid line. The major difference between the earlier L3 result and the 
revised data was the use of the average normalization from the output of two different Monte Carlos. We 
find it remarkable that the cross section factorization theorem of eq. (19), using only input from the additive 



quark model and vector meson dominance, gives a reasonable prediction of the experimental data over a 
cross section magnitude span of more than a factor of 10^ and an energy region of 3 < < 100 GeV. On the 
other hand, a literal interpretation of the experimental data at the higher energies might indicate that the 
77 cross section is rising slightly more rapidly than our prediction, a consequence perhaps of hard processes 
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not accounted for by the vector dominance model. More accurate data are required for confirmation of this 
hypothesis. 



5 Conclusions 

The available data on nn, 7p and 77 reactions lend strong experimental support to the factorization hy- 
potheses of 

• the well-known cross section factorization theorem of eq. and eq. jis|), 

• the lessor- known nuclear slope factorization theorem of eq. ( p^ ) and eq. (pO|), 

• the relatively obscure requirement of eq. ( ^ ) that p^^ = p^p — p^^, 
as well as 

• verifying the additive quark model by measuring k — 0.661 ±0.008, a result within « 1% of the value of 
2/3 expected for the quark model, using B^p measurements over a wide span of energies, 3 < -y/s < 200. 

• confirming vector dominance using dnn, o'7p and cr-y^ over an energy region of 3 < ^/s < 100 GeV and 
a cross section factor of over 10^. 

The QCD-inspired modcl[^ that we use fits the pp and pp data on total cross sections, p- values and B 
quite well and thus gives a good phenomological fit to those data. We emphasize that the conclusions 
on factorization that we presented above are rooted in the available high energy experimental data for 
nn, 7p and 77 collisions and do not depend on the details of the model used to fit nn data. 



A Appendix 

For completeness, we summarize here the formulae and parameters needed to calculate nucleon-nucleon 
scattering, taken from ref. |^ and ref. which should be consulted for more detail. 

We model the gluon-gluon contribution to the nucleon-nucleon cross section following the parton model 

(yga{s)^Cgg j Y.ggQ{s~ml)Fgg{xi,X2)dT , (23) 

where S^g = Q'Kaj./m'^ is a normalization constant, s — ts, and 

^gg = / / fgixi)fg{x2)S{T - X1X2) dxi dx2 ■ (24) 



Note that for the parameterization of the gluon structure function as Jg{x) = Ng{l ~ a:)'Ya^ , where 
Ng — ^(6 — e)(5 — e) • • • (1 — e)/5!, we can carry out explicitly the integrations, obtaining 



= CggSggiVg 2^\——i ^0 I + — ^iog(To) 1 1 (25) 

where C^g = Cgg/9, tq = rng/s, a(0) = -a(5) = -411/10, a(l) = -a(4) = -975/2, a(2) = -a(3) = -600, 
6(0) = h{h) = -9, 5(1) = &(4) = -225, and 6(2) = 6(3) = -900. The normalization constant C'gg is a fitted 
parameter and the threshold mass mo is determined by experiment. The role of mo, which is the onset of 
(Tgg(s) with s, is somewhat analagous to the role played by p™'" in the minijet models. However, numerical 
exercises show that the value of mo is not dependent on energy and that the fit is not very sensitive to the 
value. 

Also our quark-quark and quark-gluon cross sections will be parameterized following the scaling parton 
model. We approximate the quark-quark contribution with 

7i 



agg{s)^Y.gg[C + CZ7ge^J , (26) 
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where C and G'^^^g^ are parameters to be fitted. The quark-gluon interaction is approximated as 

So 



(27) 



where the normaUzation constant C]^^ and the square of the energy scale in the log term sq are parameters 
are parameters to be fitted. 

In summary, the even contribution to the cikonal is 



+S,,C^°/log-W^(6;V^7^) 

■^0 



(28) 



The total even contribution is not yet analytic. For large s, the even amplitude in eq. (g) is made analytic by 
the substitution (see the table on p. 580 of reference along with reference ||l^) s se~*'^/^. The quark 
contribution Xqq('S, b) accounts for the constant cross section and a Regge descending component (cx l/^/s), 
whereas the mixed quark-gluon term Xqg('Sj b) simulates diffraction (cx log s). The gluon-gluon term Xgg(s, 6), 
which eventually rises as a power law s^, accounts for the rising cross section and dominates at the highest 
energies. In eq. (^), the inverse sizes (in impact parameter space) /iqq and /igg are determined by experiment, 
whereas the quark-gluon inverse size is taken as ^/^qq/^gg. 



Fitted 



Fixed 



C = 5.65 ±0.14 
C^°9 = 0.103 ±0.026 
C'gg = (1.12 ±0.05) X 10- 
CSe= 25.3 ±2.0 
Codd = 7.62 ± 0.28 
So = 16.9 ±4.9 GeV2 



as = 0.5 

e = 0.05 

mo = 0.6 GeV 

Hgg = 0.89 GeV 

^igg = 0.73 GeV 

Modd = 0.53 GeV 

k = 1.349 ±0.045 



Table 1: Value of the parameters. 

In addition to the parameters in ref. the cosmic ray data fit requires the specification of another 
parameter fc, the proportionality constant between the measured mean free path (in air) and the true mean 
free path in air. The major difference between the parameters of Table |l| and the parameters of ref. Q is 
that the errors of C'gg and so are now smaller by a factor of w 2, due to the large lever arm of the hign 
energy cosmic ray points. This in turn leads to significantly smaller errors in our predictions for high energy 
cross sections, since the dominant term at high energies is (Tgg(s). 
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